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We calculate the exciton binding energy in single- walled carbon nanotubes with narrow band gaps, 
accounting for the quasi-relativistic dispersion of electrons and holes. Exact analytical solutions of 
the quantum relativistic two-body problem are obtain for several limiting cases. We show that 
the binding energy scales with the band gap, and conclude on the basis of the data available for 
semiconductor nanotubes that there is no transition to an excitonic insulator in quasi- metallic 
, nanotubes and that their THz applications are feasible. 

QJ , I. INTRODUCTION 

In the 1940s it was predicted by P. R. Wallacei that due to its honeycomb lattice, graphene, a single monolayer of 
^ ' carbon, should exhibit unusual semimetallic behavior. The gap between the valence and conductance band is exactly 
zero, and the low-energy excitations are massless chiral Dirac fermionSf^^^ On the other hand, the electronic band 
structure of other carbon-based materials shows substantial differences from the band structure of graphene. Among 
them are fuUerenes, which can be considered as zero-dimensional carbon molecules with discrete energy spectra, and 
carbon nanotubes, which are obtained by rolling graphene along a given direction and reconnecting the carbon bonds.^ 
^ [ The energy spectrum of a single- walled carbon nanotube is determined by the way it is rolled^ and is characterized 
by two integers (n,m) (0 < m < n) denoting the relative position, = nai + ma.2, of the pair of atoms on a 
' graphene strip which coincide when the strip is rolled into a tube (ai and a.2 are the unit vectors of the hexagonal 
^ . lattice). For most combinations of n and m the energy spectrum of the nanotube is characterized by the gap, whose 
' value is comparable to those in semiconductor materials. However, for a third of n and m combinations, namely when 
T— I , n — m = 3p, p = 0,1, 2, ... the value of the gap is drastically reduced and lies in terahertz frequency range. Moreover, 
for m — n the gap vanishes in zero magnetic field, and opens only after the application of a magnetic field parallel to 
lO the nanotube axis.^^- 

^sj ' Optical properties of carbon nanotubes have been investigated by many authors;^^— It was shown that the excitonic 
T— I [ effect plays an important role and that the properties of the excitonic resonance can be modulated by applying external 
fieldsJ^^— However, excitons were theoretically studied mostly for semiconductor carbon nanotubes with sufficiently 
large gapsJ^"— For metallic nanotubes mostly the excitons associated with high branches of the nanotube spectrum 
J> , separated by an energy of about 2 eV were considered J^i^ On the other hand, the analysis of the long-wavelength 
properties of narrow-band nanotubes is also of high interest, as there is a growing number of proposals using carbon 
. , nanotubes of this type for THz applications, including several schemes put forward by the authors of the present 
' worki^"— The only work on excitons in narrow-gap nanotubes, in which the stability of excitons in metallic carbon 
" " ' nanotubes subjected to Aharonov-Bohm magnetic flux was considered is Ref. [l^, where it was claimed that even in 
this case the exciton binding energy never exceeds the gap. The main conclusions of Ref. [13 are based on numerical 
calculations in the k • p scheme. In our present work we give further consideration to this interesting case, using a 
noticeably different semi-analytical approach based on a Dirac-like matrix equation. 

At first glance, narrow-gap carbon nanotubes fit ideally the excitonic insulator picture^^ predicted for the case when 
the exciton binding energy exceeds the band gap. However, the energy spectrum of a narrow-gap carbon nanotube is 
quasi-relativistic with the effective mass of the charge carriers proportional to the gap leading to a drastic reduction 
of the exciton binding energy with reducing the gap. In this paper we consider excitons formed by relativistic quasi- 
one-dimensional electrons and holes interacting via a model yet realistic potential. We show that the binding energy 
scales with the band gap and conclude on the basis of the data available for semiconductor nanotubes that there is no 
transition to an excitonic insulator in quasi-metallic carbon nanotubes and that their THz applications are feasible. 
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II. FORMALISM AND INTERACTION POTENTIAL 



The Hamiltonian for a single free electron in a narrow-gap carbon nanotube can be written (see Appendix \^ as 

where q — is the operator of the wavevector along the nanotube axis and we use the basis IV'a) , l^s) with indices 

A, B corresponding to the carbon atoms of two different sublattices in the honeycomb lattice. Here v-p is the Fermi 
velocity in graphene, connected to the tight-binding matrix element of electron hopping \t\ « 3 eV and the graphene 

lattice constant a by hvY = -^|t|a, where a — |ai| = |a2| = 2.46 A.^ For the (n,n) armchair nanotube the value 
of the band gap 2fi,WF| A| is determined by the external magnetic field, A = ^^-2^ sin , where $ is the magnetic 

s2 f Since for 



flux through the nanotube cross section, $o = ch/e is the flux quantum and h = 



I - i cos- 



experimentally accessible magnetic fields ^/^o <C 1, in our further consideration we set 6 = 1. A similar Hamiltonian 
can be written for a narrow-gap carbon nanotube with a gap opened by curvatureS^ or for certain types of graphene 
nanoribbons21. The diagonalization of equation ^ gives a quasi-relativistic dispersion, E = ±hvF\/ + q^. To go 
from the case of the electrons to the case of the holes, \'\j)f) — )■ \ipi) one should use the substitution q^—q, A^-— A. 
For a pair of interacting electron and hole the total Hamiltonian can be written in the form of a 4 x 4 matrix, and 
the stationary Schrodinger equation for determining the binding energy written in the basis = \1pf)\^p''') reads 
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where the indices e and h correspond to the electrons and holes, and q^. ^ = 
band-filling effects this Hamiltonian yields four energy eigenvalues corresponding 
particles: E = hvF{±y/ A^ -I- q'^ ± y/ A"^ + qf^); only the solution with positive sig: 
system containing a single electron and a single hole. As the potential of electron- 
only on the distance between the electron and hole, it is convenient to use the 
coordinates, X = {x^ + Xh)/2,x = x^ — Xh, q^. = K/2 + k, qu = K/2 — k, and 



as ^ij{X, x) = e 



j{x), which permits the substitution of the operator K by 



meaning of the wavevector of the exciton as a whole. Considering the case of K 
the equation for the wavefunction of relative motion reads 



. In the absence of interaction and 
to a pair of non-interacting quasi- 
ns should be taken if we consider a 
hole interaction V{xe — Xh) depends 
center of mass and relative motion 
represent the exciton wavefunctions 
the number K , having the physical 
corresponding to a static exciton, 
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where s = E/hvp, V{x) — V{x) /hvp and k = —id/dx. Eq. Q represents a system of first order differential equations, 
which can be reduced (see Appendix [B|) to a single second order equation for (f)AA'- 



^AA 



dx^ 



1 dV{x) 
V{x) dx 
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Before solving Eq. we need to specify the interaction potential. It needs to possess the following properties: first, 
it should remain finite as a; — >■ and for small x scale as V{x) « — e^/ (^eV x'^ + d^) , where e is the effective dielectric 
constant and d is the short-range cut-off parameter, which is of the order of the nanotube diameter; second, for large 
X it should decay exponentially due to the effects of screening necessarily present in nanotubes of the metallic typei^ 
- it is apparent that the size of the exciton should not exceed the mean separation between quasi-particles. The 
convenient choice of the potential is 



Vix) 



cosh(/3x) 



(5) 
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FIG. 1: The dependence of the exciton binding energy £i,//3 on the interaction strength: (a) - for a semiconductor nanotube 
with A//3 = 1; (b) - for a quasi-metalhc nanotube with A//3 — 0.01. The different lines correspond to different excitonic states. 
The crosses indicate the exact analytic results for zero-energy states (et = £g)- 



where a — / {ehvpd) and the effective length, 1//3, defines the spatial extension of the interaction. The additional 
advantage of the potential given by E g. (EI is that it allows one to obtain some analytical results in the context of 
graphene physics as was shown in Ref. [23- 



III. RESULTS AND DISCUSSION 



The results of the numerical solution are shown in Fig. [T] and Fig. [2] Figures [Ija) and [Hb) show the dependence 
of the binding energy Eh calculated as an absolute value of the difference between the eigenenergy of Eq. (jU and the 
energy of the pair of non-interacting electron and hole, £(, = |e — 2A|, measured in the units of /3 versus the effective 
strength of the interaction a//3 for two different values of A//3 — 1 and A//3 — 0.01 corresponding to the cases of 
the semiconductor and narrow-gap quasi-metallic nanotube respectively. In both cases, for small values of a//3 there 
is only one bound s-type state, whose binding energy increases with the increase of the interaction strength, until 
it reaches the value Sf, — Eg = 2A after which it goes to the continuum of states with negative energies and thus 
unbinds. On the other hand, the increase of a/f3 leads to the appearance of higher order solutions, corresponding to 
p, d, etc. excitons. This is an adequate picture for a single impurity-like exciton. However, if the interaction between 
an electron and a hole created across the gap separating the ground and excited states of the system exceeds the band 
gap in a many-body system, the ground state should be redefined and the many-body effects should govern the value 
of the gap. To find out whether the gap in a narrow-gap nanotube is mostly governed by single-electron effects such 
as an external magnetic field and curvature or by many-body eff'ects, the value of the effective interaction strength 
a/P should be estimated from the known data for excitons in semiconductor nanotubes. 

Note also, that by reducing Eq. (j?]) for e = to the hypergeometric equation (see Appendix |B]) one can find the 
analytical expression for the values of the parameter a/ 13 when the exciton binding energy is exactly equal to the 
bandgap: 




(6) 



a 

where n = 0,1,2,... . These exact values are shown by the crosses in Fig. [1] 
It is convenient to express the dimensionless interaction strength a//3 as 

a _ _ c 1 _ 300 1 

^ ~ edhvF/S ~ ^hcTpd 137 el3d' ^' 

When the interaction strength is written in the above form, one can immediately see the direct relation of the strongly- 
bound exciton problem to the long-standing problem of a supercritical charge (the nuclear charge with Z > 137) and 
atomic collapse in relativistic quantum mechanics^^. This problem has recently been revisited and reformulated for 
graphcnc'^'' where the effective impurity charge is increased by a factor of c/vf « 300, which is also present in the 
numerator of Eq. ([7]). It follows from Eq. ^ that for an exciton in a narrow-gap carbon nanotube (A ^ /3) the 
effective supercritical charge corresponding to Eb > Eg can be achieved for a//3 > 2. In principle, parameters e and /3 
can be controlled by submerging nanotubes in a solution or changing the number of electrons and holes by injection, 
optical excitation or varying temperature. However, one should not expect the value of l/{e(3d) to be higher for a 
quasi-metallic carbon nanotube than for a semiconductor nanotube of similar diameter. 
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FIG. 2: The dependence of the ratio of the exciton binding energy to the nanotube band gap on the interaction strength for a 
quasi-metallic carbon nanotube with A//3 = 0.01 (upper curve) and a semiconductor nanotube with A//3 = 1. Only s-states 
are shown. The crosses correspond to the exact analytic results when the exciton binding energy is equal to the band gap. 

Fig. [2] shows the dependence of the ratio of the binding energy of the Is exciton to the value of the gap on 
the interaction strength for the cases of the semiconductor and narrow gap nanotubes. For very small interaction 
strengths the ratio £hl^g is a universal quadratic function oiaj fi (see Appendix [C|) . The two curves remain practically 
indistinguishable from one another with increasing a//? up to aj fi 0.6. A further increase of the interaction strength 
or range leads to the breakdown of universality for finite range potentials. 

It should be noted, that for a zero-range attractive potential between the electron and hole"^^ this universality 
holds for the whole range of interaction strength (see Appendix |D]) . However, for the smooth interaction potential 
considered here the universality breaks for sufficiently large interaction strength. The difference between the case of 
semiconductor and quasi- metallic nanotubes is not as large as one could expected. Namely, the value of a//3 for which 
the Is-exciton binding energy is exactly equal to the band gap is about 3.2 for the case of a semiconductor nanotube 
with A//3 = 1 and about 2 for a quasi-metallic nanotube with A//3 = 0.01. Also, the variational calculations^^ of the 
binding energy in semiconductor nanotubes supported by the experimental data gives a value for the exciton binding 
energy of approximately 30 percent of the bandgap which according to Fig. [5] corresponds to the values of a//3 « 1. 
Using the experimental data for the exciton binding energy for the excited subbands of metallic tubes^'^ provides 
even smaller values for q;//3. These estimations ensure that for a quasi-metallic nanotube the Is exciton bound 
state lies within the gap. The hypothetical situation a/ fi > 2 corresponds to the case of the so-called "excitonic 
insulator" where the account of many body effects becomes crucially important and goes beyond the scope of our 
consideration, which is essentially a two-particle approach. Certain many-body aspects of the physics of excitons in 
narrow gap nanotubes and their relation to the Luttinger liquid model have been recently studied using a numerical 
renormalization groupi^ 

Fig. [3] shows the square modulus of the excitonic wavefunction in real space for a quasi-metallic nanotube. One sees 
that for the Is state the density in the center of mass has a local minimum, which differs strikingly from the result 
obtained earlier for semiconductor nanotubes in the effective-mass approximation^i^ in which the probability density 
in the exciton center of mass has a maximum. This is related to the complex matrix structure of the Hamiltonian 
([T]) resulting in the multi-component structure of the eigenfunctions. A similar dip in the ground-state density was 
previously reported for graphene-based waveguides;^ 

It should be noted that taking into account the valley and spin quantum numbers increases the number of different 
types of excitons associated with a given carbon nanotube spectrum branch to sixteen.^ - Their consideration, however, 
can be done along the lines described above, the only difference being the modified interaction strength a//3 for each 
type of exciton. The exciton with the highest binding energy in a semiconductor nanotube is known to be optically 
inactive (dark). The difference between the energy of the dark and bright excitons is proportional to the exciton 
binding energy and exceeds k^T at room temperature causing a significant suppression in the optical emission from 
semiconductor nanotubes. As it is shown above, for narrow-gap nanotubes the binding energy is drastically reduced 
and there should be no noticeable difference in the population of dark and bright excitonic states at any experimentally 
attainable temperatures. At room temperature, all dark and bright excitons in narrow-gap nanotubes should be fully 
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FIG. 3: (Color on-line) The density of the Is-exciton for a (10,10) carbon nanotube with a magnetic field induced gap of 
10 meV (2.5 THz) corresponding to a magnetic field of 15 T along the nanotube axis. The density represents the probability 
of finding the electron and hole comprising the exciton at the indicated relative separation. Red and blue colors correspond to 
the highest and lowest values of density, respectively. 

ionized and the direct inter-band transitions^i^— govern the emission in the terahertz range. 

IV. CONCLUSIONS 

In conclusion, we considered the formation of the exciton in narrow gap carbon nanotubes characterized by the 
quasi-relativistic spectrum of free particles. We show that the exciton binding energy scales with the band gap 
and vanishes with decreasing the gap even for strong electron-hole attraction. Therefore, excitonic effects including 
strongly-bound dark excitons, which explain the poor electroluminescent properties of semiconducting nanotubes, 
should not dominate for narrow-gap carbon nanotubes. This opens the possibility of using quasi-metallic carbon 
nanotubes for various terahertz applications. 
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Appendix A: Armchair carbon nanotube in a magnetic field 

Graphene's effective matrix Hamiltonian following the notations of Ref. |^ is given by 

H = t 



fk 
ft 



and the corresponding eigenvalues are 



£; = ±|t|VI/fel', 



where 



fk = exp (^-^k.^ -I- 2exp (^i-^^j=k^ cos • (Al) 



For an (n, n) armchair carbon nanotube, kx is quantized in the following manner 

27r I 



where I is an integer. Defining kx as the projection of the wavevector along the nanotube axis Eq. (jAl[) can be 
expressed as 

fk = exp z — - + 2exp -i-- cos — - . (A2) 



3 n J \ 3 n J \2 
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In the presence of a magnetic field, I l + T (here T = $/$o) and the armchair carbon nanotube energy spectrum 
becomes 



E^±\t\ ^1 + 4cos (-^) cos + 4cos2 (^^) . (A3) 

We are interested solely in the lowest conduction and highest valence subbands, which correspond to / = n. In this 
instance the subbands energy spectra are 



£; = ±|i|^/l-4cos(^^:^jcos(^^^ +4cos2 (^^), (A4) 

where the + (— ) sign denotes the lowest conduction (highest valence) subband. The introduction of a magnetic field 
along the nanotube axis shifts the minimum in the spectrum away from kx = 27r/3a, the electrons "acquire mass" 
and a band gap is opened. The new minimum, A:™'", is obtained by differentiating Eq. (jA4p and equating it to zero, 
this yields: 

Icosfn-) =cos(fc™-J). (A5) 



Since we are interested in particle behavior in the vicinity of fc™'", it is natural to re-express the electron energy 
spectrum in terms of qr, defined as the momentum measured relative to fc™'" i.e. kx = fc™'" + qt, thus Eq. (jA2l) can 
be expressed by 



fk = exp 



Using the identity Eq. (|A5p fk becomes: 



exp I ITT — ) + 2 cos 



[krr + qr) - 



fk = exp 



1 



T 



2J1 



1 



■ faqT\ 

sm 

V 2 / 



exp ITT 



Expanding Eq. (|A6[) in terms of qx and retaining first order terms only yields 



fk = exp 



1 + 



J" 



— 1 / 1 — T COS^ I TT 



J" 



aqr 



The effective matrix Hamiltonian can therefore be written as: 



Ho = t 



-iO 



(-*sin(^f ) - .^l-icos2 {iT^)aq^ 



''' ^sin(^f)- Jl-icos2 {■K^)aq 



(A6) 



(A7) 



(A8) 



where 9 = f (l + ^) and q 



' dy' 



This 2x2 Hamiltonian acts on a two-component Dirac wavefunction with 



components xi and X2 associated with the A and B sublattices of graphene respectively. By changing the basis 
wavefunctions from xi and X2 to ■i/'A = — Xi and V'b — ~ exp(— i0)x2, and changing the variable y — >■ —x, so that the 
x-coordinate is now along the nanotube axis, the effective matrix Hamiltonian can be expressed as 



Hq = \t 








1 — i cos^ {jr^^aq — i sin (tt^) 
1 - icos2 (7r|^)ag-|-isin(7r|^) 



(A9) 



which acts on the two-component Dirac wavefunction (iPai^b)'^- For brevity let 6 = y I ^ I ^os^ (^•^) ^^'^ ^ — 
^ sin (^.F). Thus, 



bq-iA 



bq + iA 







(AlO) 



where now q 



"dx- 
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Appendix B: Derivation of the analytic solution for E — 



For the case of total momentum K = 0, corresponding to the static exciton, the multi-component wavefunction 
of relative motion satisfies the matrix equation given by Eq. ([3]) of the main text. Since (pAA — 4'bBi the system of 
equations ([3]) reduces to 



s-V{x) 



k + lA) 4>AA 



(Bl) 



'JAB 



£-V{x) 



k-iA] (j)AA; 



BA+ [k + iA] (pAB = 



£~V{x) 



')AA- 



Substituting Eq. (jBll) and Eq. (IB2[) into Eq. (IB3|) yields the following second order differential equation: 

2 



1 dV{x) d4> 



AA 



dx^ 



£ - V{x) dx dx 



+ 



e-V{x) 



^AA = 0, 



which coincides with Eq. Q of the main text. 
Let us now consider the case of 



V{x) 



cosh (/3x) 

and e = 0. Making the change of variable z = /3x transforms Eq. (|B4I) to 



d'^(l)AA , d(j) 



tanh (z) ■ 



AA 



--^u? tanh^ iz) + '^u? - A^ 
4 4 



OAA 



= 0, 



c?z^ ' ^ ' dz 

where oj = a/j3 and A = A//3. The change of variable x = tanh(z) allows Eq. (|B5|) to be expressed as 



where c — — and e = ^cj^ — A^. Eq. (|B6p is of a known form and has the solutions^^ 



4>AA = A, (x+if (x-ir^? 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



where Ax is a constant, the function ?7 (x + 1)] is to be found and p and q are found from the following conditions: 

(B8) 



4q(<7- 1) +2(7 + c + e = 0; 
{p-q) [2{p + q)-l]=Q. 



Eqs. (Uni mi |B8|) yield 



U'-l)S + [(2p + 2g + l)x-2(p-g)]^+ [(p + gf + c 
dx^ dx L 



?7 = 0. 



(B9) 



Performing a change of variable (x + 1) /2 ^ k reduces Eq. (IB9p to the hypergeometric equation 



d^ry 

K (k - 1) + [(fli + a2 + 1) K - as] — + aia2?7 = 0, 



dr] 
d/i 
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where ai = p + g ± \/—c, 02 = p + V~-c and = 2p + i. Thus, the form of 77 is 



V^2Fi(p + q- V^,p + q + ^/^;2p+^; k] . (BIO) 



Hence the wavefunction 4>aa is given by 



0^^-Ai(l + xf (l-x)SFi (p + 9-|,P + 9+|;2p+i; i±^') . (Bll) 



From Eq. (jB8[) g is found to be 



and p can take the values 

p = q (B13) 

or 

p = i - g. (B14) 
Let us first consider the case of p = q. In this instance Eq. (|Blip becomes 



.AA = [1 - x'] ' 2F1 ( 2g - |, 2g + |; 2g + 1; ) . (B15) 



For the function c/jaa to vanish as z — )■ 00 we require 



1 ± + 4A2 ^ 
g = \ > (B16) 

and that the hypergeometric series is terminated. This can be satisfied if we take the positive root of q and restrict 
uj such that 2q — ^ = —n, where rt is a positive integer, thus we arrive at the condition that 



= 1 + 2n + + 4A2. (B17) 



The functions (pBA and (pAB: which can be found from Eq. (IBll) and Eq. (|B2I) . must also vanish as z — > cxd. For 
simplicity let us analyze the linear combinations ^+ = 0b a + 0ab and = (t>BA — (f'AB- 

^i>^ = i-———— = -i-cosh{z)——; B18 
V [z) oz UJ oz 

4A 4A 

^ _ = -i-——-(f)AA ^ i — cosh (z)0AA- (B19) 

V (Z) UJ 

For the functions satisfying Eqs. ()B18|) and (|B19p to vanish in the limit of z — > 00, we require 



.^ '"^ >1. (B20, 

This is automatically satisfied for any finite A. Therefore, providing (j^AA vanishes as z — > 00, (/>ba and (/>ab will also 
vanish. 

Let us now consider the case oi p + q = ^ , in this instance Eq. (jBlip becomes 

c^AA = ii (1 - xr (1 + X)^-'' (1 - |, i + |; ^ - 2q; i±2^) . (B21) 

Using the identitjj^i 

zi-'^sFi (a + 1 - c, 6 + 1 - c; 2 - c; z) = 2F1 (a, b; c; z) , (B22) 
where a = 2g - |, 6 = 21? + |, c = 2(7 + i and z = (1 + x) /2, Eq. ((B2Tt becomes 

= A,[1- xr^Fi (29 - 2g + |; 2q + \; i±^) , (B23) 

where A2 is a constant. This is of the same form as Eq. (jBlsp . 
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Appendix C: Binding energy in the limit of wealc electron- hole attraction 

Let us consider the limit of very weak electron-hole attraction: |V^(a;)| ^ A for any inter-particle separation. 
Substituting e = 2 A — et into Eq. (|4]) and retaining only the first order terms in V / A and eb/A transforms Eq. (|4]) 
into 

1 (P4>AA 



A +^(^)(I>AA^ -Sb(t)AA. (CI) 

Eq. (|Cip is of the form of the non-relativistic one-dimensional Schrodinger equation, where the parameter A, plays the 
role of the effective mass. For a particle of mass m in a weak one-dimensional potential U (x), the textbook^^ treatment 

r oo 1 2 

of the potential as a perturbation yields the following result for the binding energy: ~ {m/2fi^) U{x)dx , 

which can be reformulated for Eq. (|C1|) as 

-I 2 



A 

"' = J 



V{x)dx 



For V{x) = -a/cosh(/3x), Eq. ^ yields 



Sb^^A {a/l3f. 



(C2) 



(C3) 



Thus, we have shown that for a weak hyperbolic secant potential, when a ^ A, the binding energy scales with the 
band gap and has the universal quadratic dependence on a/ (5. Notably, Eq. (jC3p is also true for a weak interaction 
potential given by a Lorentzian, V{x) — a/ (l + /3'^x'^). 

Appendix D: Exciton with attractive delta- function potential 

The problem of calculating the energy of an exciton in a carbon nanotube can be solved exactly in the case when 
interaction potential between the electron and hole is taken to be a delta- function, V{x) = —UoS(x), where the 
strength of the potential Uq is positive and can be estimated as a product of the strength of the realistic potential 
and its width. 

Let us consider an electron-hole pair in a narrow-gap ID carbon nanotube with a band gap energy of 4| A|. Measuring 
all quantities in units of hvp as above, one can write the Hamiltonian as 



H = VA2+g? + JA2 +ql + V{\xe - xh\), 



(Dl) 



where 



'}e,h 



dx, 



e,h 



(D2) 



Here we have retained only the eigenstates of the general graphene-type Hamiltonian with positive energies. The 
Schrodinger equation for the problem we consider thus reads 



VA^ + lf + JA2 +ql + Vi\Xe - Xh 



^{Xe,Xh) = e'^{Xe,Xh). 



(D3) 



Introducing new variables corresponding to the center of mass and relative motion in the manner explained in the text 
of the article and putting the wavevector of the center of mass motion equal to zero, one gets the following expression 



2a A2 



d2 



Vix) iPix), 



(D4) 



where "0(2;) is a wavefunction of relative motion. The above expression is a ID Schrodinger equation for a particle with 
a complicated dispersion placed in an external potential V{x), which can be solved for the case when V{x) = —UqS{x). 
Indeed, in this case the solution of the Schrodinger equation corresponding to a bound state reads 



iPix) = y4e-''(')^ a; > 0, 



(D5) 
(D6) 
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where e = 2-\/A^ + k^. It is well known that the zero-range delta potential is equivalent to the introduction of the 
specific boundary condition for the derivative of the solution at a; = (the function itself should be continuous, 
?/'(+0) — 0)). To obtain a condition for the derivative, let us integrate the Schrodinger equation 



ipix) = UoS{x)ip{x), 



(D7) 



where Hq (^) ~ 2y + across the interval [—0; +0]; this yields 



dx 



Ho[-r] i^ix)dx\,^+o - Hoi — ] i>{x)dxU=-o = Uoi>iO) 



dx 



(D8) 



Using this relation and Eqs. (jD5[) . (|D6p and (ID4p one gets the following expression for determining the parameter k 

Hoi-K) + HoiK)-2Ho{0) = -KUo. (D9) 



Thus, 



which yields 



4 



8f7oA 
16 + 



where the energy is given by 



(DIO) 



(Dll) 



(D12) 



Note, that this result is valid only for the case of Uq < 4, as for Uq > 4, the right hand side of Eq. (|D10|) is negative. 
The binding energy of the exciton can be determined as 



Eb = 2A-e = 4AC/oV(16 + f/o'), 



(D13) 



which tends to the non-relativistic result of = [/q A/4 for Uq ^ 1; c.f. this result with Eq. (IC2|) . Note, that in both 
relativistic and non-relativistic cases the binding energy is proportional to the gap. 
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